Isospin mixing in nuclei within the nuclear density functional theory 



O 
O 

(N 



o 

=5 



> 
(N 
00 



o 

On 
O 



X 



W. Satula, 1 J. Dobaczewski, 1 ' 2 W. Nazarewicz, 3,4 ' 1 and M. Rafalski 1 

1 Institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, 00-681 Warsaw, Poland 
2 Department of Physics, P.O. Box 35 (YFL), FI-40014 University of Jyvaskyla, Finland 
3 Department of Physics & Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
* Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831 

(Dated: March 6, 2009) 

We present the self-consistent, non-perturbative analysis of isospin mixing using the nuclear den- 
sity functional approach and the rediagonalization of the Coulomb interaction in the good-isospin 
basis. The largest isospin-breaking effects are predicted for N—Z nuclei and they quickly fall with 
the neutron excess. The unphysical isospin violation on the mean-field level, caused by the neutron 
excess, is eliminated by the proposed method. We find a significant dependence of the magnitude 
of isospin breaking on the parametrization of the nuclear interaction term. A rough correlation 
has been found between the isospin mixing parameter and the difference of proton and neutron 
rms radii. The theoretical framework described in this study is well suited to describe a variety of 
phenomena associated with isospin violation in nuclei, in particular the isospin symmetry-breaking 
corrections to superallowed Fermi beta decays. 
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The isospin symmetry, introduced by Heisenberg [l[ 
and Wigner is largely preserved by strong interac- 
tions; a small violation of isospin on the hadronic level is 
due to the difference in the masses of the up and down 
quarks In atomic nuclei, the main source of isospin 
breaking is the electromagnetic interaction 0, [B[ . Since 
the isovector and isotensor parts of electromagnetic force 
are much weaker than the strong interaction between nu- 
cleons, many effects associated with isospin breaking in 
nuclei be can treated in a perturbative way. With this 
caveat, the formalism of isotopic spin is a very powerful 
concept in nuclear structure and reactions @, 0] , where 
many spectacular examples of isospin symmetry can be 
found. 

The main effect of Coulomb force in nuclei is to exert 
a long-range overall polarization effect on nuclear states 
whose detailed structure is dictated by the short-ranged 
strong force. The net effect of such a polarization is 
a result of two competing trends: the nuclear force is 
strongly attractive in the isoscalar neutron-proton chan- 
nel, while the Coulomb force acts against this attraction 
by making neutron and proton states different. In order 
to explain this interplay, self-consistent feedback between 
strong and electromagnetic fields must be considered to 
best locate the point of the nuclear equilibrium. 

An excellent example of this interplay is the system- 
atic behavior of nuclear binding energies: with increas- 
ing mass number, the stability line bends away from the 
N=Z line towards the neutron-rich nuclei. The effect of 
electromagnetic force on nuclear binding is clearly non- 
perturbative. Even in medium-mass nuclei, which are of 
principal interest in this study, energy balance between 
strong and Coulomb forces is not tremendously favorable, 
e.g., 342 MeV versus 72MeV in 40 Ca. The situation be- 
comes dramatic in superheavy nuclei and in the neutron 



star crust, where not only the binding but also spectra are 
strongly impacted by the Coulomb frustration effects re- 
sulting from a self-consistent, non-perturbative feedback 
between strong and electromagnetic parts of the nuclear 
Hamiltonian [1,0]. 

The strong motivator for studies of isospin break- 
ing is nuclear beta decay. The new data in superal- 
lowed + -^0 + nuclear beta decays [l(J, including a num- 
ber of high-precision Penning-trap measurements, re- 
quire improved calculations of isospin-breaking correc- 
tions [111, |l2]. As far as nuclear spectroscopy is con- 
cerned, there has been an increased interest in isospin- 
related phenomena in recent years 0]. For instance, stud- 
ies of excited states of proton-rich nuclei with N < Z re- 
sulted in significantly improved information on Coulomb 
energy differences 13j. In some cases, observed Coulomb 
shifts turned out to be surprisingly large [lj|, thus fuel- 
ing speculations of significant nuclear charge-symmetry- 
breaking forces. 

A precise description of Coulomb effects in nuclei con- 
stitutes a notoriously difficult computational challenge. 
In the shell-model approach to the isospin mixing (see, 
e.g., Refs. fliL flijl). the effective shell-model Hamilto- 
nian including the Coulomb interaction is diagonalized 
in a proton-neutron basis to account for non-perturbative 
effects. The overall strength of the isospin-breaking in- 
teractions is usually renormalized by reproducing the rms 
proton point radii obtained from spherical Hartree-Fock 
(HF) calculations [1611 or by fitting the experimental iso- 



baric mass shifts 




To take into account the cou- 



pling to the giant monopole resonance that appreciably 
influences the radial mismatch between the proton and 
neutron wave functions || , single-particle wave functions 
can be taken from HF calculations. More precise treat- 
ments require determining the effective Coulomb inter- 
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action in the large space, which is possible in the no-core 
shell model. Such calculations have been carried out for 
10 C [HI in the space allowing all 8 Ml excitations relative 
to the unperturbed ground state. Currently, however, ab- 
initio approaches to superallowed Fermi transitions do 
not go beyond 10 C which marks the state of the art. 

In heavier nuclei, especially those involving many nu- 
cleons outside closed shells, the isospin mixing can be 
well described by the mean-field (MF) or energy-density- 
functional (EDF) methods where the Coulomb 
force amounts to making the neutron and proton single- 
particle orbitals different, and the long-range polarization 
effects (e.g., those related to the isoscalar and isovector 
monopole resonance) are fully taken into account. 

The fact that the MF methods allow for precise treat- 
ment of long-range operators is, in fact, essential for the 
physics of i sosp in mix ing . However, it was very early 
realized 2^, 21 , 2^, 23 , 24 1 that these nice physical prop- 
erties of the MF methods are accompanied by unwanted 
spurious effects related to the fact that the neutron and 
proton single-particle states in TV ^ Z nuclei are different 
even without Coulomb interaction included. Indeed, the 
presence of the neutron or proton excess automatically 
yields isovector mean fields, i.e., different HF potentials 
for protons and neutrons. This unwelcome feature has 
hampered MF calculations of the isospin mixing beyond 
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FIG. 1: Isospin-mixing parameter otc for the even-even Ca 
isotopes determined before (BR) and after (AR) rediagonal- 
ization in the good- isospin basis \T,T Z ). The basis states 
were generated by means of self-consistent calculations with- 
out (eMF=0; upper panel) and with (eMF=e; lower panel) the 
Coulomb term. The inset shows the AR results plotted in the 
logarithmic scale. 



where n enumerates the cigenstates \n, T z ) of H . The 
value of n — 1, corresponds to the isospin-mixed ground 
state (g.s.). In the following, the g.s. isospin-mixing pa- 
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and energy E n= i^ z ob- 



the N = Z systems (see, e.g., Ref. 25]). To overcome tained after rediagonalization (AR) are distinguished 



this difficulty in the present study, we employ the mean- from the quantities a c = 1 - \b\r z \,T z \ and E\ T ^ Tz 



field methods in the framework of Refs. [22], |23( , which 
is entirely free of the spurious isospin mixing. Thereby, 
for the first time, we determine the isospin mixing within 
the context of modern EDF methods. 

We begin by noting that the self-consistent MF state 
|MF) can be expanded in good- isospin basis \T,T Z ): 



|MF) = 



T>\T 2 



b T ,T z \T,T z 



E 

T>\T Z \ 



J T,T Z | 



= 1, 



(1) 



where T and T z — (N — Z)/2 are the total isospin and 
its third component, respectively. The basic assumption 
behind our approach is that the states \T,T Z ) capture the 
right balance between strong and Coulomb interactions; 
i.e., they contain self-consistent polarization effects to all 
orders. Below we shall validate this assumption by vary- 
ing the MF charge eMF > which defines the strength of 
the Coulomb interaction at the MF level, that is, when 
solving the self-consistent HF equations. On the other 
hand, the mixing coefficients br,T z ar e not reliably deter- 
mined by the MF method, because they are affected by 
the neutron-excess-induced spurious isospin mixing. 

To assess the true isospin mixing, the total Hamilto- 
nian H (strong interaction plus the Coulomb interaction 
with the physical charge e) is rediagonalized in the space 
spanned by the good-isospin wave functions: 



\n,T z ) 



T T 



\T,T Z ), 



(2) 



T>\T Z 



(T=\T Z \,T Z \H\T=\T Z \,T Z ), obtained before rediagonal- 
ization (BR; isospin projection after variation). The AR 
results are superior to the BR results as they are based 
on the double variational principle. 

Our self-consistent calculations have been carried out 
by using the SLy4 EDF parameterization (26| and the HF 
solver HFODD [27j that allows for arbitrary spatial de- 
formations of intrinsic states. Both direct and exchange 
Coulomb terms are calculated exactly. We adopted the 
standard technique of isospin projection [28| . Details per- 
taining to our method can be found in Ref. [2§] , together 
with numerical tests. 

To illustrate the effect of the spurious isospin mixing, 
in Fig. Q] we show the BR and AR results for the even-even 
Ca isotopes. Without Coulomb interaction (eMF = 0), 
there is no isospin mixing in the N—Z nucleus 40 Ca, but 
the neutron-excess-induced mixing appears in all systems 
with N =/= Z. The BR spurious mixing is quite large, 
ac~0-2-0-4%. With the the standard Coulomb interac- 
tion (eMF = e), the BR isospin mixing increases to about 
0.2-0.7%. 

The AR results are entirely different. In 40 Ca, with 
eMF = ewe obtain the isospin mixing of 0.9%, which is 
about 50% larger than the BR value. A similar increase is 
predicted for other N=Z systems (see the upper panel of 
Fig.[2) . This result nicely illustrates the non-perturbative 
character of the Coulomb polarization when it comes to 
the isospin mixing. The impact of the isospin mixing on 
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FIG. 2: Upper panel: isospin-mixing parameter ac in even- 
even N—Z nuclei calculated in BR and AR variants. Lower 
panel: total Et=o,t z =o (BR) and E n= \ t T z =o (AR) energies 
relative to the MF (or HF) energy Bmf- 
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FIG. 4: The isospin-mixing parameter Aac (upper panel) 
and total energy AE (lower panel) in the AR method relative 
to the values obtained with the full Coulomb term, plotted 
for 40 Ca as a function of eMF- 
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FIG. 3: The AR isospin-mixing parameter ac calculated for 
even-even nuclei with 8<Z<50 and ~2<T Z <5. 



the g.s. structure of N=Z nuclei also shows up for the to- 
tal binding energy (see the lower panel of Fig. [2j ■ Differ- 
ences between the BR and AR energies rapidly increase 
with mass number, to attain about 2MeV in A=10Q. In- 
terestingly, the AR values are amazingly close to the HF 
energies Emf , up to 90 keV. This is a typical effect of the 
variational method: the minimum of energy is reason- 
ably reproduced even if the trial wave function is rather 
incorrect in its detailed structure. 

The AR isospin mixing is rapidly quenched with | N — 
Z\. Indeed, ac in a T z =l nucleus 42 Ca drops to 0.2%, 
and then decreases exponentially to about 0.01% in 60 Ca. 
As seen in Fig. this behavior holds for all isotopes. It 
is interesting to see in Fig. Q] that the AR results obtained 
for eMF = and eMF = e are quite similar beyond 44 Ca. 
This indicates that the correct good-isospin basis \T,T Z ) 
is generated from isospin-broken HF states even without 
the Coulomb interaction included on a MF level. 

This fact is further corroborated by the AR results, 
shown in Fig. Q] relative to those obtained with the full 
Coulomb term, as a function of eMF- For 0.2 e < eMF < 
e, the isospin mixing of 0.9% obtained with eMF = e 



does not vary by more than 0.01%. At eMF = 0.2 e, the 
amplitude of the \T = 2,T 2 = 0) component in the MF 
wave function becomes too small to be included in the AR 
calculation; hence, the isospin mixing jumps by 0.06%. 
This is so, because at this small value of eMF, the AT = 2 
coupling of the Coulomb force becomes ineffective and 
the numerical accuracy cannot be controlled. The |T = 
l,T z — 0) component becomes very small only in the 
very close neighborhood of eMF — 0, at which point the 
isospin mixing disappears altogether. 

The lower panel of Fig. [4] shows the total AR energy 
as a function of eMF- Here, we can understand the role 
of eMF as a variational parameter that can be used to 
optimize the good-isospin basis \T,T Z ). It is gratifying 
to see that the minimum of energy is obtained almost 
exactly at the physical value of eMF = e. Namely, the 
optimal wave functions \T,T Z ) are generated by taking 
at the MF level the full Coulomb interaction having the 
physical charge. However, it is to be noted that the en- 
ergy differences in Fig. 2] are quite small, of the order of 
a few hundred keV. Moreover, as discussed above, the 
isospin mixing is almost insensitive to such a refinement 
of \T,T Z ). This result supports our initial assumption: 
the good-isospin states \T, T z ) are fairly robust to the 
variations of the isospin-breaking interaction, i.e., they 
well capture self-consistent polarization effects. 

Although our results give first reliable estimates of the 
isospin mixing within extended MF theory, the final val- 
ues of ac are still quite uncertain, which is due to an 
imperfect determination of the nuclear EDF. This is il- 
lustrated in Fig. [5] which shows the isospin mixing calcu- 
lated in BR and AR variants for a heavy N—Z nucleus 
100 Sn for a wide selection of the Skyrme EDF parame- 
terizations [l8j ]. We note that ac does depend on the 
nuclear effective interaction: the difference between ex- 
treme AR values obtained for SkO and SkP is as large 
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FIG. 5: Isospin-mixing parameter ac in 100 Sn BR (upper 
panel) and AR (lower panel) for various Skyrme EDF pa- 
rameterizations as a function of the difference between the 
MF proton and neutron rms radii. Straight lines, represent- 
ing linear fits, are drawn to guide the eye. See Ref. [lij for 
details of Skyrme functionals used. 



as 1.5%, which is about 30% of the value of the isospin 
mixing in 100 Sn. 

In trying to pin down those features of the EDF that 
would be responsible for differences in ap, we have at- 
tempted to find correlations between isospin mixing and 
various EDF characteristics fl9| . We conclude that no 
clear correlations exist between ac and those EDF pa- 
rameters that are related to nuclear-matter properties. 
In particular, this is true for the nuclear-matter symme- 
try energy, the prime suspect to influence the properties 
of the isovector channel. We did find a very clear correla- 
tion of the BR values of ac with the differences between 
the MF proton and neutron rms radii (see Fig. [5]). This is 
not surprising, as the monopole polarization does impact 
the proton and neutron radii, and their difference. How- 
ever, after the rediagonalization, the values of ac show 
a much weaker correlation. Clearly, the precise values of 
the isospin mixing parameter depend on fine details of 
the nuclear EDF. 

In conclusion, we performed the self-consistent analy- 
sis of isospin mixing within the extended mean-field ap- 
proach. Our method is non-perturbative; it fully takes 
into account long-range polarization effects associated 
with the Coulomb force and neutron excess. The nuclear 
Hamiltonian, including the full Coulomb interaction, is 
diagonalized in a good-isospin basis obtained by isospin 
projection from self-consistent HF states. Not surpris- 
ingly, the largest isospin-breaking effects have been pre- 
dicted for N=Z nuclei, where the effects due to the neu- 
tron (proton) excess are smallest and the Coulomb force 
dominates the picture. 

The unphysical isospin violation caused by the neu- 
tron excess is significant on the MF level: the largest 
effect is predicted in \N — Z\=2 nuclei. However, the re- 



diagonalization procedure eliminates the spurious isospin 
mixing almost completely. While the optimal many-body 
solutions are obtained by using the double variational ap- 
proach, we have demonstrated that one obtains a reason- 
able good-isospin basis by broadly varying the strength 
of the Coulomb interaction of the EDF. Finally, we in- 
vestigated the dependence of isospin mixing on the self- 
consistent feedback between the nuclear and Coulomb 
terms. We found an appreciable dependence of ac on 
the parametrization of the nuclear functional and found 
a rough correlation between the isospin mixing and the 
difference between proton and neutron rms radii. The 
microscopic approach described in this study will be ap- 
plied to isospin symmetry-breaking corrections to super- 
allowed Fermi beta decays, Coulomb energy differences, 
and properties of analogue states. 
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